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Abstract 

Experiments  with  a  laminar  flow-based  membrane-less  fuel  cell  (LFMFC)  have  been  conducted  by  many  scientists.  Choban  et  al.  reported  that 
the  cell’s  performance  is  cathode  limited.  Accordingly,  we  have  established  half-cell  models  in  our  paper  to  study  two  types  of  redox  reactions 
occurring  at  the  cathode  of  the  LFMFC  without  considering  the  fuel  reaction.  Our  two-dimensional  models  are  solved  by  using  the  spectral 
method  where  the  eigenvalues  are  obtained  by  employing  the  Galerkin  method.  The  similarity  transform  is  applied  to  separate  the  concentrations 
of  the  oxidant  and  the  intermediate  product  from  their  coupled  boundary  conditions.  As  shown  in  our  results,  the  limiting  average  current  density 
increases  with  the  stoichiometry  coefficient  of  electrons  in  the  case  of  no  intermediate  product,  yet  the  maximum  electric  power  is  independent 
of  this  coefficient.  Given  the  concentrations  of  the  oxidant  and  the  intermediate  product  at  the  inlet  end  of  the  cell,  we  have  obtained  a  condition 
increasing  the  current  density.  However,  we  also  found  the  principle  of  generating  a  great  deal  of  electric  power  by  increasing  the  concentration  of 
the  intermediate  product  at  the  inlet  end  of  the  cell. 
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1.  Introduction 

The  possibility  of  generating  electrical  energy  by  continu¬ 
ally  feeding  electrochemically  active  materials  to  a  suitable  cell 
has  been  investigated  since  1839  by  using  platinum  electrodes 
immersed  in  aqueous  acid  electrolytes.  However,  important 
progress  in  fuel  cell  technology  dates  from  the  early  1950s  due 
to  the  requirements  of  the  US  Space  program  which  has  stimu¬ 
lated  many  researchers  to  seek  electrical  power  that  is  clean,  safe, 
efficient,  economically  feasible,  and  capable  of  quick  start-up. 

Many  researchers  have  developed  various  fuel  cells  [1]  for 
solving  the  two  conventional  problems  in  fuel  cells:  (1)  slow 
electrochemical  reaction  rates  result  in  low  efficiency  and  (2) 
hydrogen  has  not  become  a  readily  available  fuel  source.  On  the 
basis  of  the  electrolyte  employed  in  the  cell,  there  are  five  main 
types  of  fuel  cells:  proton  exchange  membrane  (PEM),  alkaline, 
phosphoric  acid,  molten  carbonate,  and  solid  oxide  fuel  cell.  In 
a  proton  exchange  membrane  fuel  cell  (PEMFC),  for  preventing 
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the  mixing  of  fuels  and  oxidants,  two  electrode  plates  are  bound 
to  either  side  of  the  PEM  that  is  used  for  transferring  the  protons 
from  the  anode  to  the  cathode.  Excellent  reviews  of  the  PEMFC 
up  to  mid  1990s  were  presented  by  Prater  [2]  and  Gottesfeld 
[3].  The  computational  works  on  the  PEMFC  can  be  referred  to 
the  studies  of  Bernardi  and  Verbrugge  [4,5],  Springer  et  al.  [6,7], 
Fuller  and  Newman  [8],  Nguyen  and  White  [9],  Gurauet  al.  [10], 
Yi  and  Nguyen  [11,12],  Um  et  al.  [13],  and  Kulikovsky  [14], 
However,  the  PEM  has  two  disadvantages:  (1)  it  must  be 
continuously  hydrated  to  provide  good  transport  of  protons  and 
(2)  the  fuel  can  diffuse  through  the  PEM  to  the  cathode,  caus¬ 
ing  the  problem  of  a  mixed  potential.  To  avoid  these  issues, 
Ferrigno  et  al.  [15]  and  Choban  et  al.  [16],  respectively,  devel¬ 
oped  the  laminar  flow-based  membrane-less  fuel  cell  (LFMFC) 
in  which  the  fuel  and  oxidant  aqueous  solutions  are  pressure- 
driven  from  separate  sources  into  a  common  channel.  In  this 
channel,  these  two  streams  remain  laminar  parallel  flows  with  a 
liquid-liquid  interface.  As  a  result,  the  first  disadvantage  of  using 
the  PEMFC  is  naturally  eliminated  for  the  LFMFC  because  the 
liquid-liquid  interface  is  essentially  a  non-drying  membrane. 
Moreover,  the  problem  of  mixed  potential  can  be  effectively  pre¬ 
vented  because  the  LFMFC  is  normally  operated  at  the  Peclet 
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Nomenclature 

Cint  non-dimensional  concentration  of  the  intermedi¬ 

ate  product 

Cox  non-dimensional  oxidant  concentration 

Cox  and  Cint  the  inlet-end  concentrations  of  the  oxidant 

and  the  intermediate  product  in  the  case  of  Eqs. 
(47a)  and  (47b) 

Ci  inlet-end  concentration  of  the  oxidant  in  the  case 
of  Eq.  (17) 

CBL  concentration  boundary  layer 

iav  average  current  density  generated  at  the  cathode 

electrode  by  the  reaction  of  Eq.  (17) 
t'av-i  and  iav- 2  average  current  densities  generated  at  the 
cathode  electrode  by  the  reactions  of  Eqs.  (47a) 
and  (47b),  respectively 

t'av-net  and  Pnet  net  average  current  density  and  net  average 
electric  power  generated  at  the  cathode  electrode 
through  the  simultaneous  reactions  of  Eqs.  (47a) 
and  (47b) 

/rint  and  k!ln  t  the  reaction  rate  constants  for  the  intermediate 
product  in  the  anodic  process  of  Eq.  (47a)  and  the 
cathodic  one  of  Eq.  (47b),  respectively. 
kox  and  k^d  the  reaction  rate  constants  for  the  oxidant 
(Ox)  and  the  final  product  (Rd) 

LFMFC  laminar  flow-based  membraneless  fuel  cells 
n  the  stoichiometry  constant  for  the  electron 

P  average  electric  power  generated  at  the  cathode 

electrode  by  the  reaction  of  Eq.  (17),  defined  as  r] 
time  /av 

Pe  Peclet  number 

PEM  proton  exchange  membrane 

PEMFC  proton  exchange  membrane  fuel  cells 

Greek  letters 

ri  difference  between  the  non-dimensional  electric 

potentials  of  the  cathode  electrode  surface  and  the 
adjacent  electrolyte  solution 
t]eq  the  potential  rj  at  the  equilibrium  state 


number  above  10,000  [16]  for  optimal  performance  and  for 
reducing  the  crossover  diffusions  of  fuel  and  oxidant.  Recently, 
the  principle  of  laminar  parallel  flows  has  been  successfully 
applied  to  other  microfluidic  systems  [17-20].  To  date,  a  surge 
of  studies  has  been  reported  on  micro  power  sources  for  the 
increasing  demands  of  small  and  often  portable  devices  capa¬ 
ble  of  operating  in  longer  periods  without  recharging;  e.g.,  cell 
phones,  laptop  computers. 

Ismagilov  et  al.  [21]  investigated  the  crossover  diffusion  of 
electrolytes  in  parallel  flows.  They  found  the  scaling  laws  for 
the  crossover  diffusion:  the  thickness  of  the  crossover  diffusion 
zone  is  scaled  as  the  one-third  power  of  (DZ/U)  near  the  bound¬ 
ary  and  as  the  one-half  power  in  the  middle  of  the  channel, 
where  D  is  the  diffusivity,  Z  the  distance  along  the  channel,  and 
U  is  the  averaged  velocity.  Put  succinctly,  the  thickness  of  the 
crossover  diffusion  is  determined  only  by  the  Peclet  number.  In 


summary,  these  scaling  laws  provide  a  basic  understanding  of 
the  crossover  diffusion,  which  is  required  to  predict  the  reso¬ 
lutions  of  patterning  and  the  fabrication  with  laminar  flow  [22] 
and  is  also  necessary  for  the  designs  of  diffusive  T  sensors  [23] 
and  other  components  in  micro  analytical  systems  [24] . 

Anisin  et  al.  studied  two-dimensional  models  for  the  laminar 
single  gas  layer  confined  between  a  rigid  plate  and  an  electrode 
[25].  The  electric  current,  as  shown  in  their  results,  increases 
with  the  gas  velocity. 

Bazylak  et  al.  [26]  used  a  three-dimensional  fully  numerical 
model  to  calculate  the  convective  mass  transport  in  the  LFMFC. 
Employing  the  reaction  kinetics,  Bazylak  et  al.  derived  the 
boundary  conditions  for  the  fuel  and  the  oxidant.  Their  results 
show  that  the  fuel  utilization  can  be  increased  from  8%  to 
23%  but  only  3%  of  mixing  observed  at  the  outlet  is  increased 
if  the  inlet-end  velocity  is  decreased  from  0.1  to  0.02  ms-1. 
However,  based  on  the  simulation  results,  their  proposed 
tapered-electrode  design  for  the  LFMFC  enables  the  fuel 
utilization  to  reach  more  than  50%. 

The  Butler- Volmer  equation  [27]  has  been  applied  to  calcu¬ 
late  the  current  density  generated  by  a  single  reaction.  However, 
there  are  possibilities  for  two  reactions  to  occur  simultaneously 
and  these  reactions  are  linked  by  a  common  reactant.  For  exam¬ 
ple,  the  single  reaction  O2  +  4H+  +  4e_  ^  2H2O  occurring  at 
the  cathode  electrode  can  be  regarded  as  the  sum  of  two  simul¬ 
taneous  reactions  [28]: 

02  +  2H+  +  2e“  <=►  H202  (I) 

and 

H2O2  +  2H+  +  2e“  ^  2H20  (II) 

As  seen,  reactions  (I)  and  (II)  are  not  independent  but  linked 
by  the  intermediate  product,  hydrogen  peroxide.  Experiments 
show  that  the  hydrogen  peroxide  is  a  relatively  stable  and 
detectable  intermediate  product  because  reaction  (I)  is  gener¬ 
ally  inherently  faster  than  (II).  However,  the  more  hydrogen 
peroxide  arising,  the  less  electric  current  generated  by  reaction 
(I),  because  increasing  the  chemical  potential  of  the  hydrogen 
peroxide  can  activate  the  anodic  process  of  reaction  (I).  On 
the  contrary,  because  some  of  the  hydrogen  peroxide  that 
accumulates  at  the  cathode  electrode  will  diffuse  away,  as  a 
result,  the  hydrogen  peroxide  in  reaction  (II)  is  reduced  and 
hence  less  electric  current  is  generated  by  this  reaction.  If  such 
linked  simultaneous  reactions  occur  at  the  cathode  electrode, 
the  net  current  density  generated  at  the  cathode  will  be  treated 
as  the  sum  of  those  currents  generated,  respectively,  from  the 
simultaneous  reactions  [28-30]. 

In  this  paper  we  use  theoretical  approaches  to  studying  the 
single  and  linked  simultaneous  reactions  occurring  at  the  cath¬ 
ode  electrode  of  the  LFMFC.  The  reason  why  we  focus  only 
on  the  cathode-side  reactions  is  based  on  a  report  [16]  which 
mentioned  that  the  entire  current  density  cannot  increase  signif¬ 
icantly  with  the  fuel  concentration  but  can  be  limited  obviously 
by  the  low  oxidant  concentration.  The  experiments  [16]  showed 
that  the  cell’s  performance  was  cathode  limited.  Accordingly,  in 
this  paper  we  are  interested  in  developing  two-dimensional  half- 
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cell  models  for  studying  the  cases  of  single  reaction  and  linked 
simultaneous  reactions  occurring  at  the  cathode.  We  begin  with 
the  governing  equations  for  the  steady  incompressible  parallel 
flows.  We  solve  the  mass  transport  equations  by  means  of  the 
spectral  method  where  the  eigenvalues  are  numerically  obtained 
by  using  the  Galerkin  method  [31]. 

2.  Mathematical  models 

2.1.  Continuous  and  momentum  equations 


electrostatic  force  in  the  momentum  equation  by  means  of  the 
electroneutrality  assumption  (see  Chapter  11  of  Ref.  [28]). 
Using  the  characteristic  velocity,  length,  and  pressure: 


(1) 


(**,  /)  =  H*(x,  y), 


(2) 

(3) 


Fig.  1  is  a  two-dimensional  schematic  of  the  LFMFC,  where 
the  oxidant  and  fuel  solutions  are  pressure-driven  into  a  common 
channel  and  remain  in  laminar  parallel  flows  with  a  liquid-liquid 
interface.  The  catalyst  electrodes  are  installed  at  two  bound¬ 
ary  walls  y*  =  —H  and  H* .  The  oxidant  and  fuel  solutions  are 
mainly  transported,  respectively,  in  the  cathode  side  (between 
y*  =  -//  and  0)  and  in  the  anode  side  (between  y  =  0  and 
H*).  These  two  streams  are  assumed  to  be  fully  developed 
steady  incompressible  high-Peclet-number  parallel  flows.  The 
crossover  diffusion  will  be  considered  in  the  mass  transport 
equation.  We  assume  constant  fluidic  properties  and  neglect  the 


in  which  v*  is  a  characteristic  kinematic  viscosity,  we  can  obtain 
the  non-dimensional  continuity  and  momentum  equations  for 
the  parallel  flows  in  the  common  channel: 


^=0 

3.x 


(4) 


0  = 


3  Pi  d2uj 

dx  ‘  dy 2 


0=  - 


dA 

dy 


(5) 

(6) 


Fig.  1.  Two-dimensional  schematic  of  the  LFMFC  in  which  the  oxidant  and  fuel  streams  are  pressure-driven  into  a  common  channel  and  then  remain  in  parallel 
flows  with  a  liquid-liquid  interface. 
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where  v,-  is  defined  as  v*/v*.  The  subscript  of  i  =  I  represents 
the  cathode-side  stream  that  is  transported  between  y  =  —  1  and 
0.  The  anode-side  stream  is  represented  by  ;  =  2  and  transported 
between  y  =  0  and  1 .  The  solution  to  the  velocity  u,  is  a  parabolic 
profile: 

Ui  =  - Ui[y 2  +  Ajy  +  Bt]  (7) 


Ui  = 


-dPt/dx 

2  Vj 


(8) 


We  use  four  boundary  conditions  to  solve  the  undetermined 
coefficients  in  Eq.  (7)  for  i  =  1  and  2:  they  are,  respectively,  con¬ 
tinuous  velocity  and  shear  stress  at  y  =  0  plus  no-slip  conditions 
at  y  =  ±l.  Moreover,  the  x-direction  pressure  gradients  in  the 
cathode-  and  anode-side  streams  must  be  equal  according  to  Eq. 
(6)  and  the  dynamic  boundary  condition.  As  a  result,  we  obtain: 


a  Pi  _  dP2  = 

dx  dx  P 

(9) 

8p 

Ui  =  -f- 
2  vt 

(10) 

4  M(  1  -  ru) 

A 1  =  - 

M  +  ru 

(11) 

A  1  “ru 

A  2  = 

-  M  +  ru 

(12) 

„  ru(l  +  M ) 

D  1  = 

M  +  ru 

(13) 

1  +  M 

B2  = 

M  +  ru 

(14) 

p % 

M  =  — 

P* 

(15) 

v*i 

ru  =  -r- 

(16) 

2.2.  Mass  transport  models 


In  this  subsection,  we  introduce  the  half-cell  models  and  theo¬ 
retical  approaches  to  studying  the  single  and  linked  simultaneous 
reactions  occurring  at  the  cathode  electrode.  As  described  above, 
the  fuel  reaction  is  not  considered  because  the  cell’s  perfor¬ 
mance  is  cathode  limited  [16].  The  stoichiometry  coefficients 
of  the  oxidant,  intermediate  product,  and  the  final  product  are 
assumed  to  be  one.  We  assume  high-Peclet-number  flows  for 
the  reason  described  in  the  Introduction.  In  particular,  in  the 
high-Peclet-number  flow,  the  v-direction  gradient  of  the  oxidant 
concentration  near  the  cathode  electrode  becomes  steeper,  so  a 
larger  amount  of  the  oxidant  can  diffuse  to  the  electrode. 


2.2.1.  Part  (A) 

Ox  +  2mH+  +  2 ne~  T=t  Rd  (17) 

Eq.  (17)  is  a  single  redox  reaction  where  the  notations  Ox  and 
Rd  represent,  respectively,  the  oxidant  and  the  final  product.  A 


proton  source,  such  as  sulfuric  acid,  can  be  added  to  the  cathode- 
side  stream  to  provide  a  supply  of  protons  for  the  consumption 
of  protons  at  the  cathode  electrode  [16,32,33].  In  fact,  the  con¬ 
centration  of  the  sulfuric  acid  is  very  large  as  compared  to  that 
of  the  oxidant;  the  ratio  between  the  concentrations  of  the  sul¬ 
furic  acid  and  the  oxygen  is  1000  as  shown  in  the  paper  [16]. 
Therefore,  in  such  a  very  acidic  solution,  the  concentration  of 
protons  can  be  treated  as  a  constant  [28].  In  this  paper  the  final 
product  is  assumed  to  be  pure  liquid  or  solid,  of  which  the  activ¬ 
ity  is  one.  Electromigration  is  not  considered  for  it  could  be 
small  as  compared  to  the  convective  transport  in  the  high-Peclet- 
number  flow.  Furthermore,  due  to  the  crossover  diffusion,  the 
y  variable  of  the  oxidant  concentration  must  range  from  —  1 
to  1,  so  the  oxidant  will  be  transported  by  a  different  velocity 
profile  as  it  diffuses  across  the  liquid-liquid  interface  (see  Eq. 
(7)).  Accordingly,  we  can  write  down  the  non-dimensional  mass 
transport  equation  and  the  boundary  conditions  for  the  oxidant  as 
follows: 


Peui 

9Cox 

9y 

9Cox 

9y 


(ICqx  _  92Cqx 
dx  9y2 

■  =  £oxC0xe-2^ 
=  0,  y  =  1 


-  kRd  e 


2  f'nr\ 


y=— 1 


(18) 

(19) 

(20) 


where 


Ui 


+  A,y  +  Bi ) 


(21) 


Pe  = 


—  v*SP 

n* 

uOx 


(22) 


the  superscript  of  star  (*)  stands  for  dimensional  terms  and  Cox 
is  the  non-dimensional  oxidant  concentration  based  on  the  char¬ 
acteristic  concentration  (one  molarity).  The  notation  Pe  is  the 
so-called  Peclet  number.  We  quantitatively  describe  the  reaction 
of  Eq.  (17)  by  using  the  form  of  the  Butler- Volmer  equation  as 
shown  in  Eq.  (19)  where  kox  and  kRd  are,  respectively,  the  non- 
dimensional  reaction  rate  constants  for  the  oxidant  and  the  final 
product.  Although  the  inequality  of  ft  +  f>  ^  1,  strictly  speaking, 
holds  for  the  processes  involving  more  than  one  electron  [34], 
we  will  assume  that  both  f>  and  /T  are  equal  to  0.5  in  this  paper. 
In  Eq.  (19)  77  is  a  non-dimensional  term  defined  by 


Frf 

~RT 


(23) 


where  F  is  the  Faraday  number  and  if  is  the  difference  between 
the  physical  electric  potentials  of  the  cathode  electrode  surface 
and  the  adjacent  electrolyte  solution.  Dqx  is  the  dimensional  dif- 
fusivity  for  the  oxidant.  Note  that  the  activity  of  the  final  product, 
Rd,  is  one.  To  solve  for  Cox ,  we  transform  Eqs.  (18)-(20)  as  fol¬ 
lows: 


Pe  Uj 


9Cox 

dx 


d2C  ox 
9y2 


(24) 
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9Cox 

dy 


=  eCox, 


y  =  —2 


dCox 

dv 


=  0, 


y  =  o 


(25) 

(26) 


Carrying  out  the  computation  of  Eq.  (38),  we  obtain  one  set  of 
{ai a2,k,  ■  ■  aN,k}  for  each  Ok  to  compose  the  corresponding 
i jfk  according  to  Eq.  (35). 

Let  Cox  at  the  inlet  end  of  the  common  channel  (.v  =  0)  be 
given  by 


where 


e  =  k0xe. 


—2  fnrj 


Cox  =  Cox  -  K~l  e2"" 


K~l 


^Rd 

kox 


y  =  y-  i 


(27) 

(28) 

(29) 

(30) 


C  Ox  — 


0  for  —  1  <  y  <  0 
C  i  for  —  2  <  y  <  1 


(42) 


where  C\  is  the  inlet-end  concentration  of  the  oxidant,  so  the 
inlet-end  condition  for  Cqx  can  be  found  from  Eq.  (28): 


C()x  = 


—  K  1  e2"''  for  —  1  <  y  <  0 

Ci  -  K~l  e2"''  for  -  2  <  y  <  1 


(43) 


Let 

N 

Cox  =  'ffuk  exp  fk(y)  (31) 

k=  1  e 

Substituting  Eq.  (31)  into  Eqs.  (24)-(26),  we  obtain 


By  means  of  the  orthogonality,  the  coefficients  or*  in  Eq.  (31) 
can  be  solved  by 


<*k  = 


r0  ~ 

J_ 2  C Ox(x  =  0 )ui\jfk  dv 
J_2Uifldy 


(44) 


d2fk 
dy 2 


+  rrkUii/k  =  0 


-rzr  =  efk,  aty=-2 
dy 


d  jfk 

dy 


=  0, 


aty  =  0 


We  express  i /r*  by 


N 

t k  =  y fn,k<Pn(y) 

n= 1 


(32) 

(33) 

(34) 


(35) 


To  match  the  boundary  conditions  Eqs.  (33)  and  (34),  we  may 
select 


4>n  (y)  =  COS 


(36) 


So  must  satisfy 

fin  tan  }in  =  2s  (37) 

The  eigenvalues  fin  can  be  solved  numerically.  By  substi¬ 
tuting  Eq.  (35)  into  Eq.  (32),  an  eigenvalue  equation  based  on 
the  Galerkin  method  is  obtained  as  follows,  from  which  we  can 
solve  for  a„jc  and  cr*: 

tHr1[R]fl*  =  oiajt  (38) 


where 

ak  —  [fll,*,  02,  ON.kV  (39) 


After  that,  Cox  can  be  obtained  from  Eq.  (31),  and  hence 
Cox  is  solved.  Now,  the  local  current  density  generated  at  the 
cathode  electrode  can  be  calculated  by  substituting  the  Cox  into 
Eq.  (19).  Let  the  current  density  be  rendered  non-dimensional  by 
using  the  characteristic  current  density  (FDqx1  molarity)///*. 
We  define  an  average  current  density  as 

L 

2n(k0xCox  -  kRd  e2^)|y=-2  dx  (45) 


where  L  is  the  non-dimensional  length  of  the  cathode  electrode. 
According  to  Eq.  (19),  we  have 


2 n  [LdC ox, 

lav  —  .  /  1 5>= — 2  cLx:. 

L  Jo  dy 

(46) 

2.2.2.  Part  (B) 

Ox  +  mH+  +  ne~  Int, 

(47a) 

Int  +  mH+  +  ne~  ^  Rd 

(47b) 

Now,  we  consider  two  simultaneous  reactions  linked  by  an 
intermediate  product  (Int)  as  shown  in  Eqs.  (47a)  and  (47b), 
and  their  sum  is  exactly  the  single  reaction  of  Eq.  (17),  for 
instance,  the  case  of  O2  +  4H+  +  4e_  —  2EhO.  In  this  paper 
we  assume  the  ratio  between  the  diffusivities  of  the  oxidant  and 
the  intermediate  product  is  roughly  equal  to  one,  because  the  dif- 
fusivity  for  most  electrolytes  in  dilute  solutions  is  in  the  order  of 
10~9  m2  s-1  [35].  In  the  following,  we  present  the  mass  transport 
equations  and  the  boundary  conditions  for  the  concentrations  of 
the  oxidant  and  the  intermediate  product  (Cox  and  Qm)  on  a 
basis  of  the  same  characteristic  scales  and  assumptions  as  those 
used  in  part  (A): 


\H]i,j  =  /  Ui4>i(pj  dy 


'—2 


[*]/,;•  = 


*,7=1 . N. 


(40)  .  d 

re  ui  — 

Cox 

_  92 

c  Ox 

dx 

Clnt 

dy2 

Qnt 

(4D  ^  =koxC0xe-^-kIntCInte^,  y  =  -2 

dy 


(48) 

(49a) 
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dC\n 

By 


=  ^OxCox  -  klntClnt  -  k'Cinl 


+  km^n\  y=-2 


(49b) 


3  C  ox 

351  Cint 


=  0, 


)>  =  0 


(50) 


where  Pe  and  5,-  are  defined  in  Eqs.  (21)  and  (22). We  can  com¬ 
bine  Eqs.  (49a)  and  (49b)  into: 


3 

Cqx 

'  Jkoxe-A"1? 

Cox 

By 

Cint 

-kQx  e~ kim  s^nr,+k[nt  er^nr> 

Cint 

0 


Note  that  each  electron  transfer  coefficient  in  Eq.  (5 1 ),  just  as 
like  as  /3  and  ft'  in  Eq.  (19),  is  assumed  to  be  0.5.  The  notations 
A'int  and  k\M  represent  the  reaction  rate  constants  of  the  inter¬ 
mediate  product  in  the  anodic  process  of  Eq.  (47a)  and  in  the 
cathodic  process  of  Eq.  (47b),  respectively.  However,  we  can 
rewrite  Eq.  (51)  as  follows: 


vectors  corresponding  to  A.  i  and  A.2;  that  is, 


'An 

A 12 

P\,m 

=  km 

P\,m 

A2i 

A22 

Pl.m 

Pl,m 

If  we  define 


Si 

=  ip]-' 

Si 

vr 

=  ip]-' 

h 

C()x 

Clnt 


0 

-kRd^nr> 


and 


h 

V 

1 

/i/(eki) 

h 

S2 

1 

filial) 

m  =  1,2 


(60) 


(61) 

(62) 


(63) 


then  we  can  rewrite  Eqs.  (48),  (50)  and  (52),  respectively,  as 
follows: 


Pe  Hi 


3 

Bx 


1 _ 

__  32 

l - 

•*— 

1 _ 

ti 

By2 

1 

to 

(64) 


3 

By 


=  0,  y  =  0 


(65) 


3 

Cox 

'  1  -Kienr> 

Cox 

3 

ti 

eX\t\ 

£\t\ 

By 

Cint 

—  G 

-1  ki  e'"'  +  k2 

Cint 

By 

h 

eXiti 

e2t2 

where 


0 

-*Rde*" 


As  seen,  both  t\  and  t2  are  not  coupled  as  Cox  and  Cint  in  Eq. 
(52)  (52).  The  conditions  of  t\  and  t2  at  the  inlet  end  of  the  common 

channel  (x  =  0)  can  be  found  through  Eqs.  (6 1 )— (63)  where  the 
inlet-end  conditions  of  Cqx  and  C\nt  are  given  as  follows: 


s  =  k0x 

(53) 

hnl 

Ki  =  r~ 

kox 

(54) 

-s2|  0 

^  1 

III 

(N 

(55) 

1 n 
d 
III 

(56) 

Define  the  matrix  in  Eq.  (52)  as  [A]: 


'An 

A 12 

1 

—k\  e",? 

A21 

A22 

-1 

/ci  e"'7  +  k2 

Using  the  similarity  transform,  we  have 
[A]  =  [P][A][P]-‘  (58) 


and 


[A] 


Xi  0 

0  7-2 


(59) 


where  An  and  A.2  are  the  distinct  eigenvalues  of  [A].  The  first  and 
second  columns  in  the  [7*]  matrix  are,  respectively,  the  eigen- 


CqxU=0 


0  for  —  1  <  y  <  0 

Cox  lor  —  2  <  y  <  1 


ClntU=0 


0  for  —  1  <  y  <  0 

Cint  for  —  2  <  y  <  1 


(67) 

(68) 


Therefore,  we  can  separately  solve  Eqs.  (64)— (66)  for  t\  and 
ti  by  means  of  the  same  method  for  solving  Eqs.  (24)-(26), 
and  hence  Cox  and  Cint  can  be  obtained  through  the  reverse 
transform  from  Eq.  (61)— (63).  Now,  we  can  calculate  the  net 
average  current  density  which  is  the  sum  of  two  average  current 
densities  generated,  respectively,  by  the  simultaneous  reactions 
of  Eqs.  (47a)  and  (47b).  They  are  expressed  as 

L 

n(koxCox  e^'"1  -  kmtCmt  e^)|y=_2  d.T  (69) 

L 

n(4tCint  e“^  -  kRd  ePnr,)\~y=-2  dx  (70) 

where  fav_i  and  iav- 2  are  non-dimensional  quantities  based  on  the 
same  characteristic  current  density  used  in  part  (A).  Thus  the  net 
average  current  density  is  defined  by 


Jav-net  —  1 av- I  +  *av-2- 


(71) 
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3.  Results  and  discussion 

In  this  section  we  present  the  modeling  results  of  both  current 
density  and  electric  power  as  a  function  of  potential  i]  for  the 
reactions  in  parts  (A)  and  (B).  As  defined  in  Eq.  (23),  t]  is  the 
difference  between  the  non-dimensional  electric  potentials  of  the 
cathode  electrode  surface  and  the  adjacent  electrolyte  solution. 
For  saving  on  computation,  the  potential  t]  is  restricted  to  the 
range  from  0  to  less  than  30.  Consequently,  the  conversion  by 
using  Eq.  (23)  shows  that  the  range  of  the  physical  potential  i] 
is  within  0-0.775  V  for  room  temperature  (RT/F  =  0.02585  V). 
Unless  otherwise  specified,  the  following  parameters  are  fixed: 
(1)  the  number  of  expansion  terms  in  Eq.  (31)  N=  150,  similarly 
for  t\  and  ti  expansions,  (2)  the  Peclet  number,  Pe=  106,  (3)  M 
and  ru  in  Eqs.  (15)  and  (16)  are  set  to  be  one  and  (4)  all  electron 
transfer  coefficients  in  Eqs.  (19),  (49a)  and  (49b)  are  specified 
as  0.5. 

3.1.  Results  for  part  (A) 


played  in  the  figure;  accordingly,  we  choose  10~8  and  10-19  as 
the  range  of  kRd  in  Fig.  2a  for  pure  graphic  consideration.  In  the 
equilibrium,  the  current  densities  generated  by  the  cathodic  and 
the  anodic  processes  of  Eq.  (17)  are  equal  such  that  the  deriva¬ 
tive  of  Cox  respective  to  y  at  the  cathode  electrode  is  equal  to 
0  according  to  Eq.  (19),  whereby  we  can  infer  Cqx|v=_ 1  =  Cj. 
Accordingly,  we  can  obtain  from  Eq.  (19)  the  potential  t]  at  the 
equilibrium  state: 

hi(£0xCi/A-Rd)  ln(Ciff) 

1„  = - Yn - =  j  =  ^0  (72) 


which  is  independent  of  kRd  as  shown  in  Fig.  2a.  Furthermore,  we 
combine  Eqs.  (30),  (45)  and  (72)  to  rewrite  the  average  current 
density  as 

L 

2nkRde2npri ^ 

Cox,v=-  1  e— 2/3/i(f7— »;eq )  _  e2/S«(/J-)Jeq) 

Cl 


Fig.  2a  shows  the  average  current  density  zav  versus  the  poten¬ 
tial  rj  for  the  reaction  of  Eq.  (17)  for  four  cases  of  kRd,  in  which 
we  specify  three  parameters:  the  ratio  kox  to  kRd  (K),  stoichiom¬ 
etry  coefficient  n,  and  the  inlet-end  concentration  of  the  oxidant 
(Ci).  As  mentioned  above,  we  set  the  range  of  ij  as  from  0  to  20 
(rf  =  0-0.517  V).  As  a  result,  considering  the  clarity  of  Fig.  2a, 
we  present  only  four  cases  of  fcRd  within  this  potential  range. 
Those  curves  for  kRd  larger  than  10-8  and  less  than  10“ 19  are  so 
close  to  the  curves  of  kRd  =  10“  8  and  10“ 19  that  they  are  not  dis- 


Apparently,  the  potential  ;;  must  be  decreased  much  further 
from  the  equilibrium-state  value  in  order  to  generate  appreciable 
current  density  in  the  case  of  smaller  kRd,  because  we  know  that 
the  smaller  the  kRd  and  kox  ,  the  less  the  exchange  current  density. 
For  instance,  to  generate  the  average  current  density  of  0.4,  the 
reduction  of  i]  is  found  to  be  1%  and  97%  of  rjeq  for  the  cases  of 
kRd  =  10-8  and  10“ 19,  respectively.  Now,  we  define  an  average 
electric  power  by  P  =  )/iav.  Fig.  2b  shows  that  the  maximum  P 
decreases  as  kRd  is  diminished. 

In  contrast  to  Fig.  2  where  kRd  is  varied  yet  n  is  fixed,  Fig.  3a 
shows  the  effects  of  n  on  the  cell’s  performance,  where  kRd  is 
fixed,  while  Ci  and  K  remain  the  same  as  in  Fig.  2.  As  seen, 
the  equilibrium-state  value  of  the  potential  i]  decreases  with 
increasing  n:  rje q  =  20,  10,  and  6.67  for  n=  1,2,  and  3,  respec¬ 
tively.  However,  this  property  can  be  checked  by  substituting 
these  n  values  into  Eq.  (72).  On  the  contrary,  the  limiting  aver¬ 
age  current  density  increases  with  n.  In  fact,  as  the  potential  i] 
is  lowered  to  a  certain  extent,  the  cathodic  process  of  Eq.  (17) 
overwhelms  the  anodic  one  such  that  the  oxidant  at  the  cath¬ 
ode  electrode  surface  (v  =  —  1 )  is  nearly  depleted.  Consequently 
the  y-direction  gradient  of  the  oxidant  concentration  at  y  =  —  I 
reaches  the  limit,  which  can  be  estimated  as  follows: 

dC Ox  I  0)x,bulk  region  —  ('Ox.y=— !  Cox,bulkregion  C ]  (HA \ 

dv  y=~  1  dc  <7::  dc 

where  dc  is  the  so-called  thickness  of  the  concentration  boundary 
layer  (CBF)  of  the  oxidant.  This  thickness  is  proportional  to  the 
reverse  of  the  square  root  of  the  Peclet  number  according  to  Eq. 
(18),  yet  independent  of  n.  In  the  situation  of  Eq.  (74),  the  y- 
direction  diffusion  of  the  oxidant  reaches  the  maximum  so  the 
average  current  density  will  reach  its  limit  as  well.  We  combine 
Eqs.  (30),  (46)  and  (74)  to  express  this  limiting  average  current 
density  as 


Fig.  2.  Average  current  density  zav  and  average  electric  power  P  vs.  the  potential 
r]  where  K=e 40  x  103,  n=  1,  and  C\  =  10-3.  The  symbol  of  rectangular  box 
stands  forthe  case  of  ^Rd  =  10-8;  black  circle:  &Rd  =  1 0- 12 ;  triangle:  k Rd  =  10-16; 
star:  &Rd  =  10~19. 


*av 


2 n  fL  dCpx 

L  Jo  dy 


ly=— 1  d-r  = 


L 

2nC\  I  dx 

L  J  dc 


(75) 


o 
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(b)  P 

Fig.  3.  Average  current  density  iav  and  average  electric  power  P  vs.  the  potential 
!)  where  K=em  x  103,  i.-Rd=  10-12,  and  Ci  =  10-3.  The  symbol  of  black  circle 
stands  for  the  case  of  n  =  1 ;  rectangular  box:  n  =  2;  black  triangle:  n  =  3. 


One  can  see  that  this  limiting  value  is  theoretically  linear 
with  n  because  clc  is  independent  of  n.  As  shown  in  Fig.  3a, 
the  limiting  /av  is  found  to  be  1.18,  2.37,  and  3.55  for  n  =  1,  2, 
and  3,  respectively.  Obviously  their  ratio  is  nearly  1:2:3,  which 
coincides  with  Eq.  (75).  Next,  we  display  in  Fig.  3b  the  results 
of  the  average  electric  power  P  calculated  for  those  three  cases 
of  n  in  Fig.  3a.  However,  the  maximum  P  is  about  the  same, 
nearly  independent  of  n.  Moreover,  the  potential  rj  at  which  the 
maximum  P  is  generated  is  found  to  be  13.4,  6.8,  and  4.5  for 
n=  1,  2,  and  3,  respectively.  The  ratio  among  these  potentials  is 
about  1:1/2: 1/3. 

3.2.  Results  for  part  (B) 


Fig.  4.  Average  current  densities  /av.  i  and  (av- 2  vs.  the  potential  r]  for  n=  1, 
ko\  =  102,  A:r(i  =  10  \  Cox  =  0.1,  and  Cint  =  0.  The  symbol  of  black  circle 
stands  for  the  case  of  kint  =  1 0  3  and  kf  =  108;  black  triangle;  k\nt  =  108  and 
=  108;  rectangular  box:  f|nt  =  1 0  8  and  k'^t  =  10”8. 


When  studying  the  simultaneous  reactions  of  Eqs.  (47a)  and 
(47b),  we  are  especially  interested  in  the  effects  of  the  interme¬ 
diate  product  on  the  cell’s  performance.  Consequently,  in  this 
subsection  our  discussion  is  focused  on  the  theoretical  under¬ 
standing  of  the  principles:  how  the  current  density  is  changed 
with  the  reaction  rate  constants  of  the  intermediate  product  (lcint 
and  Tjnt)  and  changed  with  the  concentration  of  the  intermediate 
product.  In  order  to  qualitatively  describe  the  effects  of  A'int  and 
A:Jnt,  in  Figs.  4  and  5,  we  investigate  three  sets  of  k\nt  and  k[nl  in 
which  we  specify  the  kox  and  k^,  and  the  inlet-end  concentra¬ 
tions  of  the  oxidant  and  the  intermediate  product  (Cox  and  Cint). 
Next,  in  Figs.  6  and  7,  we  choose  three  cases  of  Cint  to  study 
the  effects  of  the  concentration  of  the  intermediate  product  on 
the  cell’s  performance.  It  should  be  noticed  that  the  ranges  of 
/qnt  and  k'lnl  (both  are  from  10-8  to  108)  are  specified  for  pure 


Fig.  5.  Net  average  current  density  4v-net  and  net  average  electric  power  Pnet  vs. 
the  potential  rj  for  n  =  1,  kox  =  102,  k Rd  =  10-5,  Cox  =  0.1,  and  Cint  =  0.  The 
symbol  of  black  circle  stands  for  the  case  of  &int  =  10-8  and  &[nt  =  108;  black 
triangle:  A:int  =  108  and  &[nt  =  108;  rectangular  box:  k^t  =  10-8  and  k^  =  10-8. 
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Fig.  6.  Average  current  densities  idv. i  and  iav_ 2  vs.  the  potential  rj  for  n=l, 
kox  =  102,  =  10-5,  f|nt  =  1(T8,  k\M  =  10s,  and  Cox  =  0.1.  The  symbol  of 

rectangular  box  stands  for  the  case  of  Cjnt  =  0;  star:  Ci„t  =  0.01;  black  circle: 
O' [nt  =  0.05. 


30 


0  100  200  300  400  500  600  700 


(b)  Pnel 

Fig.  7.  Net  average  current  density  4v-net  and  net  average  electric  power  Pnet  vs. 
the  potential  ij  where  n  =  1,  kox  =  102,  km  =  10-5,  kint  -  10-8,  fj[]t  =  108,  and 
Cox  =  0.1.  The  symbol  of  rectangular  box  stands  for  the  case  of  C ox  =  0;  star: 
Cint  =  0.01;  black  circle:  Cint  =  0.05.  The  symbol  of  black  triangle  represents 
the  iav  and  P  calculated  for  the  single  reaction  of  Eq.  (17). 


computational  consideration,  not  for  any  particular  intermediate 
species,  because  the  main  purpose  in  this  subsection  is  to  qual¬ 
itatively  describe  the  effect  of  the  intermediate  product  on  the 
cell’s  performance. 

First,  we  display  in  Fig.  4a  the  results  of  calculated  iav_i  by 
using  Eq.  (69)  for  Aqnt  =  10-8  and  108,  whereas  £{nt  is  fixed  at  108. 
It  indicates  that  iav_i  becomes  larger  in  the  case  of  smaller  k\nt. 
In  contrast  to  Fig.  4a,  Fig.  4b  shows  the  /av. 2  calculated  by  using 
Eq.  (70)  for  k'lnt  —  10-8  and  108,  where  &int  is  fixed  at  10-8. 
As  seen,  the  4v-2  increases  obviously  with  k[nV  Yet,  it  must  be 
noticed  that  (av-2  is  almost  0  or  negative  in  the  case  of  A:Jnt  =  10~8 
because  this  k'lnl  is  a  great  deal  smaller  than  the  specified  value 
of  A’Rd  and  also  because  the  Cint  is  0  such  that  the  reaction  of 
Eq.  (47b)  is  dominated  by  the  anodic  process  especially  when 
the  potential  >7  is  enhanced.  Accordingly,  we  infer,  it  is  favorable 
to  increasing  the  net  average  current  density  4v-net  when  both 
koJhnx  and  £jnt/kRd  are  very  large.  To  prove  this,  we  present 
the  calculated  iav-net  versus  the  potential  i]  in  Fig.  5a  for  those 
cases  in  Fig.  4.  As  shown,  the  results  agree  with  our  inference. 
On  the  contrary,  the  case  of  larger  k\M  coupled  with  smaller  Cjm 
will  be  more  adverse  to  the  generation  of  current  density,  which 
is  not  our  interest  and  not  presented  in  Fig.  5.  In  addition,  we 
display  in  Fig.  5b  the  net  average  electric  power  Pnet  defined 
by  iav-net 7-  In  summary,  at  given  Cox  and  Cint,  simultaneously 
increasing  the  ratios  fcox/^int  and  k[nl/  &Rd  will  be  favorable  to  the 
enhancements  of  the  net  average  current  density  and  the  electric 
power. 

To  more  understand  the  effect  of  the  intermediate  product  on 
the  cell’s  performance,  in  addition  to  researching  the  effects  of 
A’lnt  and  k[nv  we  also  want  to  study  the  effect  of  the  concentration 
of  the  intermediate  product.  To  carry  out  this  work,  we  select  the 
optimal  case  in  Fig.  5  where  k\ni  =  10-8  and  k'lnt  =  108  because, 
as  described  in  the  last  paragraph,  this  case  is  favorable  to  the 
generations  of  current  density  and  electric  power.  Furthermore, 
we  set  Cint  =  0,  0.01,  and  0.05,  respectively.  First,  as  shown  in 
Fig.  6a,  the  calculated  iav_  1  by  using  Eq.  (69)  is  about  the  same 
among  these  three  cases  of  Cint  because  the  given  ratio  koJk\nl 
is  so  large  that  the  second  term  in  the  integration  of  Eq.  (69), 
as  compared  to  the  first  term,  is  almost  neglected  for  Cint  less 
than  0.05.  Moreover,  as  displayed  in  Fig.  6b,  the  ;av_ 2  increases 
remarkably  because  the  specified  value  of  kJnt/^Rd  is  so  large 
that  slightly  increasing  the  Cint  can  obviously  enlarge  the  first 
term  in  the  integration  of  Eq.  (70).  Accordingly,  we  infer,  only 
when  the  ratios  &ox/£int  and  A:[nt/^Rd  are  simultaneously  very 
large  will  the  net  average  current  density  obviously  increase  as 
Cint  is  increased.  To  prove  this,  as  defined  in  Fig.  5,  the  net 
average  current  density  and  the  net  average  electric  power  for 
those  cases  in  Fig.  6  are  calculated  and  exhibited  in  Fig.  7a  and 
b,  respectively.  As  seen,  both  ;av-net  and  Pnet  obviously  increase 
with  Cint.  In  particular,  the  maximum  Pnet  changes  from  156  to 
653  as  Cim  is  slightly  increased  from  0  to  0.05.  Checking  Fig.  6, 
one  can  see  that  the  evident  increases  in  iav-net  and  Pnet  are  due 
to  the  changes  in  iav-2-  Finally,  we  present  in  Fig.  7a  and  b  the  fav 
and  P  (black-triangle  symbol)  calculated  for  the  single  reaction 
of  Eq.  (17)  to  compare  with  the  4v-net  and  Pnet  at  the  same  values 
of  n,  k()K ,  kRd,  and  the  inlet-end  concentration  of  the  oxidant.  As 
shown  from  the  result  of  comparison,  adding  more  intermediate 
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product  into  the  linked  simultaneous  reactions  of  Eqs.  (47a)  and 
(47b)  can  help  us  obtain  higher  power  than  that  obtained  though 
the  single  reaction  of  Eq.  (17)  as  long  as  the  ratios  k(),Jk\nl  and 
ki„t/kRd  are  very  lar§e- 

4.  Conclusions 

Two-dimensional  half-cell  models  were  established  for 
studying  the  current  densities  and  electric  power  generated 
through  the  reactions  of  Eqs.  (17),  (47a)  and  (47b)  that  occur  at 
the  cathode  of  the  LFMFC.  The  parabolic  solutions  to  the  fully 
developed  steady  incompressible  high-Peclet-number  parallel 
flows  were  obtained  without  considering  the  electrostatic  forces. 
Mass  transport  by  electromigration  was  neglected  as  compared 
to  the  convective  transport  because  we  assumed  high-Peclet- 
number  flows  in  the  study.  We  focused  on  the  cathode-side  redox 
reactions  without  considering  the  fuel  reaction  because  a  report 
[16]  mentioned  that  the  cell  was  cathode  limited.  We  used  the 
Butler- Volmer  equation  to  describe  the  boundary  conditions  of 
the  concentrations  of  the  oxidant  and  the  intermediate  product  at 
the  cathode  electrode.  We  solved  the  mass  transport  equations 
by  means  of  the  spectral  method  where  the  eigenvalues  were 
numerically  obtained  through  the  Galerkin  method. 

For  the  single  reaction  of  Eq.  (17),  given  ratio  kox/kRd  and 
the  inlet-end  concentration  of  the  oxidant,  the  potential )]  at  the 
equilibrium  state  does  not  change  with  k rj  if  n  is  fixed,  but 
equilibrium-state  >)  decreases  with  increasing  n  if  k^d  is  fixed. 
The  latter  can  be  checked  by  substituting  these  n  values  into  Eq. 
(72).  In  addition,  in  the  limiting  condition,  the  average  current 
density  is  linear  with  n,  while  the  maximum  average  electric 
power  is  almost  independent  of  n.  However,  the  potential  t]  at 
which  the  maximum  average  electric  power  is  generated  in  the 
cases  of  n=  1,  2,  and  3  has  the  ratio  nearly  equal  to  1:1/2: 1/3. 

In  the  linked  simultaneous  reactions  of  Eqs.  (47a)  and  (47b), 
the  /'av_i  becomes  larger  in  the  case  of  smaller  k\nt  while  /av_ 2 
increases  with  k'lnV  At  given  Cox  and  Cint,  simultaneously 
increasing  the  ratios  fcox/^int  and  £jnt/A:Rd  is  favorable  to  the 
enhancements  of  the  net  average  current  density  and  the  net 
average  electric  power.  However,  the  /av_  1  is  about  the  same  for 
Cint  less  than  0.05  if  /.'ox/^int  is  very  large.  Nevertheless,  the  iav- 2 
can  increase  remarkably  by  slightly  increasing  the  Cint  provided 
that  the  ratio  £jnt/kRd  is  very  large.  In  summary,  given  the  same 
values  of  n,  kox,  /'ltd,  and  the  inlet-end  concentration  of  the  oxi¬ 
dant,  by  adding  a  great  deal  of  the  intermediate  product  into 
the  oxidant  solution,  one  can  obtain  more  power  through  the 
linked  simultaneous  reactions  of  Eqs.  (47a)  and  (47b)  than  that 
obtained  though  the  single  reaction  of  Eq.  (17)  as  long  as  the 
ratios  kox/kint  and  k'lnJ  k^d  are  very  large. 
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